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Abstract. Recently, some Lemaitre-Tolman-Bondi metrics have been considered as models 
alternative to the dark energy within the Friedmann-Lemaitre-Robertson- Walker universes. 
One of these models has been confronted with cosmic observations, this confrontation seeming 
to rule it out. On the other hand, we show that the intrinsic energy of this model does not 
physically vanish. Thus we conclude that this model is actually a non creatable one. This 
suggests from the very beginning that the model is not a good model to represent our Universe, 
which would be in according with the negative result of this confrontation with observations. 
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1 Introduction 

Out of the particular case where A' ^ (see next section), the Lemaitre-Tolman-Bondi 
(LTB) metric family is the most general family of spherically symmetric metrics in General 
Relativity corresponding to a pressure-less matter source [1]. 

Some of these metrics have been used, in a cosmological context, to describe large inho- 
mogenous structures and the anisotropies they produce on the cosmic background radiation 
temperature [2-4] and, more recently, as models alternative to the dark energy Friedmann- 
Lemaitre- Robertson- Walker (FLRW) universes [5-1 1] . The models have been confronted with 
cosmological observations, and the final conclusion of one of these papers [11] is that the con- 
frontation has become sufficiently constraining to "rule out the whole class of adiabatic LTB 
models", by testing a particular LTB model hereafter called the constrained Garcfa-Bellido- 
Haugb0lle (CGBH) model. 

On the other hand, from the beginning of the last seventies, people have speculated 
on the possibility that the Universe could have raised from a vacuum quantum fluctuation 
[12, 13], an idea further developed by Vilenkin [14]. If this had been the case, we could expect 
that the energy of our Universe, P°, and also the corresponding linear 3-momentum, P l , and 
angular 4-momentum, J a ^ , vanish ... = 1,2,3 and a,/3, ... = 0,1,2,3). In [15], we have 
called such a universe with vanishing 4-momenta a creatable one. But it is well known that 
the energy and momenta of an space-time in General Relativity are dramatically dependent 
on the coordinates used. So, what coordinates must be used in order to calculate the specific 
two 4-momenta, P a , J a ^ , that have to vanish in the case of a creatable universe? Our answer 
in [16] is that these specific coordinates have to be intrinsic ones, defined as follows: 

1) First, they are Gauss coordinates in some finite region of the considered space-time, 
covering the 3-space boundary. That is, in this region, the metric components go a take the 
values (/oo = — 1) 90i = 0. As it is well known, these are coordinates related to free falling 
observers. 

2) Second, the corresponding linear 3-momentum, P l , and angular 3-momentum, J 13 , 
vanish. 

3) Third, let it be the space-like 3-surface t = to = const., with t the time coordinate. 
We denote this 3-surface by S3. Then, asymptotically, the 3-space metric, gij, approaches fast 
enough a manifestly conformally flat metric (g^ = G 2 5ij) when we approach the boundary, 
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T*2, of S3. Of course, if the space-time is asymptotically Minkowskian, the conformal factor 
is the unity. 

Notice that in Gauss coordinates the time coordinate, t, is the proper time of those 
falling observers. Furthermore, is a universal time. This means that equal distant readings of 
this time correspond to events which are physically simultaneous (see [17], epigraph 84). This 
could be important in order to define a consistent energy, since, as it is well known and we 
are going to see, its expression (see next (1-1)) is a 3-volume integral taken in a given instant 
t. That is, the elementary partial contributions to this integral are taken for the same t value, 
and it would be a good thing that this common instant labelled physically simultaneous 
events. Therefore, it would be good to deal with a proper and universal time, i.e., with 
Gauss coordinates, when calculating what we call here the intrinsic P° energy value: by 
definition, the one calculated in intrinsic coordinates. As far as the energy, is concerned, 
the term intrinsic refers to the fact that this energy is, first, calculated for coordinates whose 
corresponding linear and angular 3-momenta, P l and J*- 2 , vanish and, second, since these 
coordinates are associated to free falling observers (local inertial ones) that, in the present 
case of LTB metrics, are further co-moving with the source pressure-less matter. Thus, these 
observers add nothing extra to the considered space-time (compare this situation with, for 
example, the static observers in a static metric that would have to be prevented of free falling 
by some virtual non gravitational action). Therefore, we choose just this coordinates to define 
the specific P° energy and momenta, the intrinsic 4-momenta, that must vanish if a given 
metric is creatable. In the present case, if some LTB models are creatable. 

These intrinsic coordinates can be proved to exist always [16] for each constant value 
t = to and, as mentioned above, the corresponding linear and angular 4-momenta P a and J a/3 
will be called intrinsic momenta. In all, we have called creatable universes the space-times 
whose intrinsic linear and angular 4-momenta, the last one irrespective of the momentum 
origin, vanish. 

We will make explicit the existence of intrinsic coordinates in the particular case of the 
LTB and FLRW models considered in the present paper. As we will see next, the closed 
and flat FLRW universes have vanishing intrinsic momenta, while the value of the intrinsic 
energy of the open FLRW universe is —00. Thus, in our parlance, the two first ones are 
creatable universes, but the last one is not. Furthermore, the perturbed flat FRW universes 
in the frame of standard inflation in neither creatable, while the perturbed closed one remains 
creatable [18]. 

Let us precise that the starting point for our definition of the intrinsic momenta of 
an space-time is the Weinberg energy-momentum complex [19]. There are many complexes 
present in the literature on the subject, but in [20] we have explained why we have selected, 
among all them, the one from Weinberg, from which some specific expressions for the energy, 
P°, the linear 3-momentum, P i , and angular 4-momentum, J a P, can be deduced. These 
expressions involve 3-space volume integrals that in the particular case of P° read 



where the gravitational constant has been taken equal to 1, g = Sijgij and the summation on 
repeated indices is performed with the Kronecker 5. This defined P° coincides (once writing it 
as a 2-surface integral on the above £2 boundary) with the well known Arnowitt-Deser-Misner 
(ADM) energy [21]. 
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Assuming that an existing universe has to be creatable is an attempt of saying something 
about our Universe all along its existence from just after the big bang. It can also supply 
us with a criteria for, from the beginning, discarding tentatively all non creatable metrics in 
General Relativity as good candidates for good cosmological models. This is just the case in 
the present paper where the CGBH model of the Universe is observationally discarded. This 
discarding could have been suggested from the very beginning because of the physically non 
creatable character of this model. We will prove this non creatable character along the present 
paper. 

2 The LTB metrics and the general expression for their energy 

The LTB metric class [1] can be written as [10]: 



satisfying the Einstein field equations, with A' = d r A and k<l. This family of metrics 
describes the spherical solutions of the Einstein field equations with a pressureless matter 
source [1]. We can see that the coordinates used are adapted to this spherical symmetry. 
Furthermore, they are Gauss coordinates, that is, for the metric components go a we have 
goo = —1, goi =0, such that the observers associated to these coordinates are free radially 
falling observers co-moving with the pressureless matter source. 

We will assume that A' exists and is different from zero everywhere, except perhaps for 
r = 0. Note that from (2.1) we recover the FLRW metrics by putting A = a(t)r, k = nr 2 , 
with a the corresponding expansion factor and k the curvature index, k = 1,0, —1. 

Because of the manifest spherical symmetry of the metric (2.1), the linear 3-momentum, 
P l , and the angular 3-momentum, J y , associated to this metric, vanish. Thus, the coordinates 
used in (2.1) will be intrinsic coordinates, as defined above, provided that the r coordinate be 
such that dl 2 becomes in the boundary of the corresponding 3-space a manifest conformally 
flat metric. Then, we are left with the question whether the energy, P°, calculated in such 
coordinates, that is, the intrinsic energy, vanishes or not in order to conclude if a particular 
LTB metric is creatable or not. 

According to some general expressions given in [22], the expression for the P° energy of 
the metric (2.1) becomes 



where we have put both the gravitational constant and the speed of light equal to 1. 

To obtain this expression one must transform the 3-space integral giving P° in (1.1) into 
a 2-surface integral on the boundary of this 3-space by applying Gauss theorem. To apply this 
theorem we need that the metric be regular enough: the 3-space derivatives of the 3-space 
metric must be continuous. 1 Thus, we will assume not only that, except for r = 0, A' in (2.1) 

1 If there are no intrinsic singularities in the integration 3-volume, there always exist coordinates in which 
the Gauss theorem can be applied, provided that we assume, as it always done, that the differentiable manifold 
of the General Relativity is C 2 -class by pieces. Nevertheless, in the present work, we only have to use intrinsic 
coordinates. Then, the C 2 -class by pieces character could not be fulfilled if we are restricted to only use these 
intrinsic coordinates. Therefore, for intrinsic coordinates, the above regularity condition in order to apply 
Gauss theorem should be verified each case for the particular metric used. 



l-Jfc(r) ' v ' 

with da 2 = d6 2 + sin 2 9d(j) 2 and where A and k are functions of the corresponding arguments 





(2.2) 
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exist everywhere and is different from zero, but further than A" is continuous too, and that 
there is no intrinsic singularity at r = 0. 

As remarked above, since we want this energy P° to be an intrinsic one, we must use in 
(2.2) an r coordinate such that dl 2 in (2.1) be asymptotically conformally flat in a manifest 
way. 

In the next section, we calculate P° for two different versions of a particular LTB metric, 
called above the CGBH model [10, 11], which are regular enough in the above precise sense. 



3 The CGBH model 



In references [10, 11], the authors explore the possibility that we live close to the center 
of a large void, i.e. close to the center of a suitable LTB model, as an alternative to the 
prevailing interpretation of the Universe acceleration in terms of a ACDM model with a 
dominant dark energy component. They confront this possibility with a series of cosmological 
observations through two versions, the flat and the open one, of the CGBH model cited 
above, the first (second) version becoming asymptotically a flat (an open) FLRW without 
cosmological constant. The CGBH model, in its two versions, is ruled out as a result of the 
confrontation. 

Because of this asymptotic character, comparing the LTB metric (2.1) with its flat FLRW 
universe limit for r — > oo, having in mind the expression ds 2 = — dt 2 + a(t) 2 (dr 2 + r 2 da 2 ) for 
this flat FLRW metric, we easily obtain 

A(t, r -)• oo) ~ a(t)r, lim Mr) = 0. (3.1) 

Notice that, because of this asymptotic character, the coordinates in (2.1) are asymptotic 
conformally flat coordinates for dl 2 . Consequently, they are intrinsic coordinates (they fulfill 
the three conditions 1 — 3 defining the notion of intrinsic coordinates in the Introduction) and 
the corresponding P° energy we are going to calculate will be the intrinsic energy. 

Because of (3.1), the partial contribution to P°, in (2.2), coming from the term containing 
the k function is 

1 rkA' 2 a 2 rk a 2 

o *!P i — u = ~ir I? 1 i — u = ( r/c )' ( 3 - 2 ) 

Z r— >oc 1 — K Z I — >oo 1 — K 2 r— >oo 

whose actual value depends on how fast k vanishes when r — > oo. Then, in accordance with 
[10], let us define the function Qjif(r), that generalizes the matter cosmic parameter, 17, of 
the FLRW cosmology, by writing k{r) like 

k = A 2 (n M -l), (3.3) 

with Aq = dtA{t = to, r) and to the cosmic present time. 

The value of this function for the flat version of the CGBH model is [10] 

\l M = il out + 1J — — T7 j-r-T , (3.4) 

1 + tanh(ro/ Ar) 

with Q ou t = 1 (£lout < 1) for the open case), where f2j n , ro and Ar are parameters to be fitted 
by cosmological observations. In particular, ro characterizes the void size, near whose center 
we are assumed to be placed, and Ar the transition to uniformity. 



-A- 



From (3.4), we easily obtain 



n M (r» Ar) - 1 ~ Ae" r/Ar , (3.5) 

with 

2 e ro/Ar 

= 1 , , r /A (n»n-l)- (3.6) 
1 + tann ro / Ar 

Then, from (3.3) and the first equation of (3.1) we have 

k(r » Ar) ~ Xa 2 r 2 e- r/Ar . (3.7) 

According to (3.2), this asymptotic behavior of k means that the contribution to P° in 
(2.2) from the term involving k vanishes. Thus, in order to calculate P° for the flat CGBH 
model, we are left with the remaining contribution from the term dealing with the function 
A. To calculate this contribution we have to make explicit how fast A approaches a(t)r (see 
the first equation of (3.1)) when r — > oo. The easiest way to make this is to use one of the 
Einstein field equations for the LTB metrics (2.1) (see [10]): 

A 2 + 2AA + k = 0, (3.8) 

2A/A + All A' = -4vrp A/ . (3.9) 

More precisely, we will consider Eq. (3.8) for r — > oo, where as we have just seen k 
behaves like k ~ r 2 e~ r / Ar , jointly with the parametric A value for the flat CGBH model 
[10, 11], i.e., 

A(t,r) = Qm (cosh ri- l)r, (3.10) 
2(1 - Um) 

2(1 - tt M r /z 

with rj the parameter, and 

Ho(r) = §Ur (l " Om)- 1 ^ " n o M u/2 sinh-^i^)^], ( 3 .12) 
2 (1-SZm) 7 I'M 

with the inserted factor 3/2 allowing us to obtain the Hubble constant Hq, for Hq(t —> oo). 
Finally, Q,m is given by (3.4) with O, ou t = 1- 

Then, in accordance to (3.1), let us write the asymptotic form of A like 

A(t,r -too) = a(t)r[l + e(t,r)], (3.13) 

where e is a function of t and r such that lim^oo ~ = for any t value and for any p > 1. This 
asymptotic expression for A would guarantee the vanishing of the corresponding contribution 
to the intrinsic value of P°. In the Appendix we show that A has actually such an asymptotic 
behavior since we obtain 

e(t,r^oo)~e- r/Ar . (3.14) 
Then, the resulting asymptotic form for A when r — > oo leads to 
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, (A-rA') 2 

lim ^ '— = (3.15) 

and so to the final vanishing of the intrinsic energy of the flat CGBH model. 

It seems, then, that the asymptotically flat CGBH cosmological model is creatable. 
However, we are going to see in the Section 5 that this vanishing of P° is a dramatically non 
stable result. We will account for this limitation by saying that the model is physically non 
creatable. But, before considering the instability question, in the next section we will obtain 
the asymptotic form of the functions A(r,t) and k(r) in (2.1), that is, for a general LTB 
metric, when the r coordinate allows us to write dl 2 asymptotically in a manifest conformal 
way. 

4 LTB metrics in asymptotic conformally flat coordinates 

Let us make a transformation of the radial coordinate r in (2.1), going to a new radial co- 
ordinate p, trough a time independent function, that is r — > p = p(r). We will chose this 
function such that the 3-space metric, dl 2 , becomes in the new radial coordinate asymptoti- 
cally conformally flat in a manifest way. That is, 

dl 2 (p -> Pb ) ~ G 2 (t,p)5 ij dp i dp j = G 2 (dp 2 + p 2 da 2 ), (4.1) 

where pi is such that 5ijp l pi = p 2 and where p = p b = constant is the equation of the radial 
boundary of the space-time given by the metric ds 2 in (2.1), this equation having be assumed 
time independent. We can expect such a radial coordinate to exist since in [15, 16] it has been 
proved on general grounds for every constant time t = to. This existence is obvious in the 
particular case of the CGBH model since, as it is mentioned at the beginning of the Section 
3, this model approaches asymptotically an open or flat FLRW universe. 
Then, by simply comparing (4.1) with (2.1), we obtain 

A'(t, p -> p b f/[l - k(p -> p b )] ~ G 2 (t, p -> p b ), (4.2) 
p 2 b G 2 (t,p^p b )~A 2 (t,p^p b ). (4.3) 
We will write these two equations more compactly by putting 

A' 2 /(l - k) ~ G 2 , A 2 ~ p 2 G 2 . (4.4) 
These asymptotic equations give trivially: 

k ~ 1 -p 2 A' 2 /A 2 , (4.5) 

which, since k only depends on p, means that the function A(t, p) factorizes asymptotically 
in its two arguments, that is to say, in an obvious notation: 

A(t,p)~a(t)f(p). (4.6) 
Having in mind this factorization, Eqs. (3.8) and (3.9) become asymptotically 

a 2 + 2ad = p 2 f' 2 / / 4 — 1/ f 2 (4.7) 

and 
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4 

a/a = --Trp M (t), (4. 



respectively, with pu{t) = p M (t,p^ Pb) 
From (4.7) we obtain 



a 2 + 2ad = C, p 2 f 2 /f - 1/f = C, (4.9) 

with C an arbitrary constant. 

Eq. (4.8) and the first one equation of (4.9) are equivalent to the two dynamical cosmic 
equations for the expansion factor, a(t), of a FLRW with —C curvature, i.e., equivalent to 
(4.8) jointly with 



d\2 _ 8irpM C 
a) 3 a 



) =^ + - 2 . (4.10) 



As far as the second equation of (4.9) is concerned, its general solution is a family of 
solutions depending on an arbitrary constant, laying aside the constant C. In all, when the 
parameter C vanishes we obtain the values for a(t) and f(p) corresponding to a flat FLRW 
universe, that is, a(t) = const, x t 2//3 , f(p) = p, while for C/Oa particular solution is 



/ = P (i - \c P 2 y\ (4.ii) 



giving for G the expression: 



G 2 {t,p) = a 2 {t){l-\C P 2 )- 2 , (4.12) 
i.e., the FLRW universes with — C curvature, in conformal flat 3-space coordinates: 

ds 2 = - d ^ + a\t) ^^f w . (4.13) 

For C > (C < 0) we have the open (closed) model and for C = the flat one. Actually, 
from (4.13) we obtain the corresponding standard form 

ds 2 = _ dt 2 + a 2 (t) ^_^_ + r 2 d(J 2^ ; (4 M) 

by making the coordinate transformation 

r = p (l - Cp 2 /A)-\ (4.15) 



Obviously, in the particular case that C > (the open case) , we could also have obtained 
the result (4.14) directly from the A(t,r) solution of the open CGBH model, that is from 
(3.10)-(3.12) and (3.4) with n out < 1. 

Finally, in the open case, normalizing C to 1, we see that, while r goes from zero to 
infinite, p goes from zero to 2, so that the above boundary value p^ becomes pb = 2. 



-7- 



5 The instability of the energy vanishing: discussion and conclusion 



Let us calculate the energy P° of any LTB metric that behaves asymptotically like an open 
(non flat) FLRW metric as we have discussed in the precedent section. To make this calcu- 
lation we will start from Eq. (1.1). Then, since we have assumed at the beginning of the 
Section 2 that our LTB metrics are regular enough, we can, using Gauss theorem, write (1.1) 
as the 2-surface integral on the corresponding boundary: 



P° = — lim f (djga - d i g)n i p 2 sm0d9d6, (5.1) 

l07T p-S-Pb J 



with n* = p l / p. 

Since now we have asymptotically dl 2 ~ G 2 (t, p)bijdx % dx^ , P° becomes 

P° = ~ lim f d i G 2 n i p 2 sm0d6d(t> = -\ lim p 2 d p G 2 , (5.2) 

07T p-+pb J 2 p->Pb 

that is, according to (4.12), with C normalized to the value (7=1, 

P° = -2limp 2 4- (1 ~p 2 /^r 2 = -oo, (5.3) 
p^2 dp 

Therefore, P°, whose value for the flat CGBH model (Section 3) was zero, jumps to a 
minus infinite value, P° = — oo, when an elementary shift of the C constant from its original 
value is performed, or what is the same, when a shift from Vt ou t = 1 (see (3.4)) to the new 
value Vt ou t < 1 is performed as close to 1 as we want. 

Thus, the vanishing of P° for this flat CGBH model is unstable and so this vanishing 
could be considered un-physical. Then, though in our parlance the flat model is strictly 
speaking a creatable one, we can consider and denote it, jointly with the open CGBH model, 
as physically non creatable universes. All this could be seen in accordance with the fact 
that this model is ruled out by its confrontation with cosmological observations [10, 11]. In 
conclusion: had this confrontation not still taken place, in view of this physical non creatable 
character, we could have predicted a subsequent negative result, not as a true prediction but 
as some plausible suggestion. 

Furthermore, assuming a negative value for the constant C in (4.13) and normalizing to 
C = —1 we obtain the closed FLRW metric in 3-space conformal flat coordinates: 

ds 2 = -dt 2 + a 2 (t) ^0J^ . (5.4) 

The boundary limit pb is now pi, = oo. Then, using the expression (5.1), we straightfor- 
wardly obtain P° = 0, that is, the LTB metric approaching asymptotically the closed FLRW 
model is creatable. Thus, we could ask if this "closed" LTB model would not be better entitled 
than the CGBH one to be observationally tested as a model of the void universes considered 
in the present paper. Actually, aside with other models, the "closed" LTB model has been 
tested in [9], without fully conclusive results. But, obviously, the creatable character of this 
model, by itself, could not avoid that finally the model were ruled out by observations. In 
fact, as stated in [9], "in practice we could always approximate the correct answer by setting 
f2fc to a small nonzero value". 

Finally, let us come back to Eq. (5.2). Notice that putting 
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lim p 2 d p G 2 = 
p^pb 



(5.5) 



fully characterizes the subfamily of creatable LTB metrics. For instance, if pi, = oo, G 2 ~ -p, 
p > 1, for p — > pb, jointly with the particular case in which G does not depend on p, describe 
all the corresponding LTB creatable universes. 



A Appendix. Calculating the e(t, r) function 

We start from (3.8), that is, A 2 + 2AA + k = 0, and from (3.7), that is, k(r » Ar) ~ 
\a 2 r 2 e- r / Ar . Then we write A as A = a(t)r[l + e(t,r)] with a(t) = (t/t ) 2/3 . After an 
elementary calculation, neglecting quadratic terms in e, we find 

-^(i 2 + 2AA) ~ (d 2 + 2aS)(l + 2e) + 2a(ae + 3de). (A.l) 
But we have a 2 + 2ad = 0. Then, Eq. (3.8) becomes for large values of r 

2a{ae + 3de) ~ -dgAe~ r / Ar , (A.2) 
that in accordance with a(t) = {t/to) 2 ^ can be written as 

(a 3 e)^~Xe- r / Ar , (A.3) 
whose general solution for large values of r (r >> Ar) is 



e(t,r) = f{r)a(t) + h(r), 



(A.4) 



with 



f(r) 



^Ae- r / Ar , 

9t 2 



h{r) an arbitrary function, and a(t) the general solution of 



(a a) = a, 



that is to say, 



« = ^t 4/3 t 2/3 + pt^r 1 + v, 



(A.5) 



(A.6) 



(A.7) 



with p and v two arbitrary constants. 

Substituting this expression of a in (A.4), we obtain, for large values of r but for any 

time, 



e(t,r) 



(A. 



Then, in order to fix both arbitrary constants, p and u, and the arbitrary function h(r), 
let us come back to the CGBH model, more precisely, to Eqs. (3.10) and (3.11). For small rj 
values, that is, for small t values, we obtain 
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A(t,r) ~ (^) 2/3 ^ 3 (r)^ 2/3 (r)rt 2 / 3 , (A.9) 

where we can substitute Ojvf by its asymptotic expression from (3.5). 

Furthermore, having in mind (3.4) with Q ut = 1 and (3.12), after an elementary calcu- 
lation, we obtain the following asymptotic value 

H(r» Ar)~H {l-l\e- r / Ar ). (A.10) 

5 

Carrying this expression to (A.9) we obtain for large values of r and small values of t 

A(t, r) ~ ar [l + -(1 - a) A e~ r/Ar ] . (A.ll) 
5 

Finally, through the relation A = ar(l + e), let us compare this approximated expression 
of A(t, r) with e(t,r), given by (A. 8). We obtain for large values of r, but for any time, 

e~ ^A(l-a)e- r / Ar , (A.12) 
5 

which is in accordance with Eq. (3.14). 
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